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Abstract 

In a trapped atomic Fermi gas, one can tune continuously via a Feshbach resonance the effective pairing interaction between 
fermionic atoms from very weak to very strong. As a consequence, the low temperature superfluidity evolves continuously from 
the BCS type in the weak interaction limit to that of Bose-Einstein condensation in the strong pairing limit, exhibiting a BCS- 
BEC crossover. In this paper, we review recent experimental progress in atomic Fermi gases which elucidates the nature of the 
r^J- superfluid phase as the interaction is continuously tuned. Of particular interest is the intermediate or crossover regime where the 
i— I ,.5-wave scattering length diverges. We will present an intuitive pairing fluctuation theory, and show that this theory is in quantitative 
agreement with existing experiments in cold atomic Fermi gases. 
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1. Introduction 

Ultracold atomic Fermi gases have been a very exciting, 
rapidly developing field, which has emerged within the past 
several years, bridging condensed matter and atomic, molecular 
and optical physics [1]. Using a Feshbach resonance in a mag- 
netic field, one can tune the effective pairing interaction strength 
between fermionic atoms from very weak to very strong yfl. 
! As the interaction strength varies, the nature of the low temper- 
ature superfluidity of these Fermi gases evolves continuously 
from the BCS type in the weak coupling limit to Bose-Einstein 
condensation (BEC) in the strong pairing limit, exhibiting a 
BCS-BEC crossover, which has been of great theoretical inter- 
est since 1960's . Of particular interest is the unitary 
regime, where the s-wave scattering length diverges. This is a 
strongly correlated regime where no small parameter is avail- 
able for perturbative expansions. It has been expected that this 
regime provides a prototype for studying both high T c super- 
conductors and strongly interacting Fermi gases which are also 
of interest to nuclear and astro-physicists. 

In this paper, we first review experimental progress in atomic 
Fermi gases, with an emphasis on recent radio frequency spec- 
troscopy measurements. Then we will present a pairing fluctua- 
tion theory compare with experiment. We show that this theory 
successfully explain experimental measurements. 

2. Experimental progress 

The first theoretical study of BCS-BEC crossover dates back 
to 1960's, although it did not get much attention until the sem- 
inal work of Leggett in 1980 on BCS-BEC crossover at zero 
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temperature J3]. The discovery of high Tc superconductivity in 
1986 gave a strong boost to the interest in BCS-BEC crossover 
It was suggested that the unusual pseudo- 
gap phenomena in the cuprate superconductors might have to 
do with BCS-BEC crossover. Experimental efforts in this area 
fell far behind, because it had been difficult to find a system 
where the attractive pairing interaction is tunable. Thanks to 
the laser cooling and trapping technique in 1990's, one is able 
to create "artificial" many-body systems of fermionic atoms in 
a laboratory. The existence of a Feshbach resonance in these 
Fermi gases makes it possible to tune the interaction strength. 

For ease of control, the Feshbach resonances for the two 
widely studied species, 6 Li and 40 K, are both very wide. The 
interaction in both cases are of the short-range, s-wave type. 
They are often taken to be a contact potential in theoretical 
treatments. 

The first experimental realization of BCS-BEC crossover 
was achieved in 2004 by Jin and coworkers, iTTTl Fl2h and almost 
the same time by the Grimm group 11311 and the Ketterle group 
lfl4ll . Due to the difficulty in tuning temperature T, the Fermi 
gases were either in the superfluid or normal state at given inter- 
action strength (or the magnetic detuning). Continuous varia- 
tion of the system as a function of temperature was first realized 
by the Thomas group II 1 511 at unitarity. In collaboration with the 
theory group at Chicago lfl6h . Thomas et al %vh observed for 
the first time continuous phase transition from the normal to su- 
perfluid state in a unitary 6 Li gas. One could argue, of course, 
that the vortex measurement of the Ketterle group provided the 
most definitive smoking gun for a superfluid state. lfl8ll 

Besides the interaction strength, another great tunability is 
population imbalance between the two fermionic species to be 
paired lfl9ll . It adds a whole new dimension to the phase dia- 
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gram and makes the physics much richer. It also generates in- 
terest [20] in possible observation of the Larkin-Ovchinnikov- 
Fulde-Ferrell (LOFF) state 112 111 . Experimental work in popula- 
tion imbalanced Fermi gases was pioneered by the Hulet group 
1 22] and the Ketterle group |23]. Experiment in the extreme 
population imbalanced limit by the Ketterle group manifested 
1 24] the importance of Hartree-like correlation effects besides 
BCS-type of pairing. 

Unlike an electron system, it has been difficult to measure the 
excitation gap in the Fermi gas superfluid. Among all experi- 
mental techniques, Radio frequency (RF) spectroscopy [ 13J is 
arguably the most direct probe. Using a tunable RF field to ex- 
cite one of the two pairing atoms from a lower hyperfine state 
(level 2) to a higher hyperfine level 3 which do not participate 
in pairing, a higher frequency will be needed if the atoms in 
level 2 are paired. Such a frequency shift (detuning) provides a 
good measure of the excitation gap. Previous measurement by 
Grimm and coworkers lfl3ll . and later repeated by the Ketterle 
group 12411 . was done in a momentum integrated fashion. At 
low T, the RF spectra displayed double-peak structure, with a 
sharp peak at zero detuning and a broad peak at positive detun- 
ing. This double-peak feature was nicely interpreted 1251 12611 
as transitions from unpaired atoms the trap edge (correspond- 
ing to the sharp peak) and from a distribution of paired atoms 
in the inner part of the trap (broad peak). However, doubt was 
cast about the origin of the two peaks as to whether they reflect 
pairing of bound state effects 12411 or simply a result of trap in- 
homogeneity 12711 . Recently, attention was also drawn to final 
state effects both theoretically I28I.I29I1 and experimentally 1 24 ] . 

A big step in the RF technique was the recent momentum- 
resolved RF spectroscopy experiment in 40 K by the Jin 
group i30ll . With momentum resolution, RF spectroscopy is 
equivalent to the angle-resolved photoemission spectroscopy 
(ARPES) for an electron system, In fact, it is cleaner than 
ARPES in that ARPES is only a two-dimensional probe, which 
is often plagued by the existence of surface states, surface con- 
taminations, work function, and the complication of energy 
dispersion in the third dimension. In comparison, of course, 
the signal-to-noise ratio in a Fermi gas experiment is much 
lower, as limited by the (low) total number of atoms in the 
gas. although the trap inhomogeneity adds complication to 
the interpretation of the spectrum. Like ARPES, momentum- 
resolved RF spectroscopy measures the fermion spectral func- 
tion, A(k, to), which is of central importance in characterizing 
the system. 



3. Theoretical Formalism 

In this section, we now present a simple pairing fluctuation 
theory, which was first developed fioll to explain the pseudo- 
gap phenomena in high Tc superconductors. Fermi gases in the 
presence of a Feshbach resonance can be effectively described 
by a two-channel model J2|]. It has now been known that the 



closed-channel fraction 13 lL 13211 is very small for both 6 Li and 
40 K, throughout the BCS-BEC crossover. Therefore, for these 
systems, a one-channel model is often used as a good approxi- 
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Figure 1: Schematic diagrams for the fermionic self-energy The dotted 

and (red) double lines represent the condensate and finite momentum pairs, 
respectively. 



mation, given by the grand canonical Hamiltonian 



+ Zq,k,k' U(k, k ')a q/2+k , T fl q /2-k,| fl q/2-k',| a q/2+k 



(i) 



where e k = h k /2m is the free atom dispersion. The difference 
between Eq. ([D and its BCS counterpart is that BCS keeps only 
the q = term in the interactions. The inclusion of finite q 
terms allows incoherent, finite momentum pairing. For clarity 
of presentation, we will take a contact potential, t/(k, k') = 1, 
and use a 4-momentum notation, K = (k, iu) n ), Q = (q, /Q/), 
Yjk - T lZklZn> an d set ?j = 1. Population imbalance can be 
described by p^ + p[. However, here we will only present the 
equations for the case of equal spin mixture. Generalization to 
population imbalance can be found in Ref. i33ll . 

We assume that (i) the fermionic self energy £ has a pairing 
origin, (ii) pairs can be either condensed or fluctuating with a fi- 
nite momentum, and (iii) condensed and noncondensed pairs do 
not mix at the level of T-matrix approximation. Figure 1 shows 
diagrammatically the contributions to the self-energy, where the 
double (red) lines indicate finite momentum pairs and the dot- 
ted line indicates the condensate. The subscripts "sc" and "pg" 
stand for superfluid condensate and pseudogap contributions, 
respectively. 

To tackle this problem, we use a Green's function method. 
We derive the equations of motion for one- and two-particle 
Green's functions G and G2, which will involve higher order, 
three particle Green's functions G3: iG = [H, G] ~ G, G2, 
iGi = [H, G2] ~ G,G2,G3. We then truncate the equations 
of motion at the level of G3, factorize G3 into a sum of prod- 
ucts of G and G2, and treat G and G2 on equal footing. For 
G2, we focus on the particle-particle channel, neglecting the 
particle-hole channel which normally only provides a chemi- 
cal potential shift. We emphasize that it is the particle-particle 
channel that gives rise to superfluidity. After some lengthy but 
straightforward derivation, we obtain the self energy: 
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MX) = -KGo(-K), 
L pg (K) = J]t pg (Q)G (Q-K), 



where 



t Pg (Q) = J] T 



u 



+ Ux(Q) 



(2) 
(3) 
(4) 



(5) 



is the (pseudogap) T-matrix, and x(Q) = TjK Go(Q - K)G(K) 
is the pair susceptibility. Here Go is the bare Green's function. 
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Figure 2: Typical behavior of (a) the temperature dependence of the gaps in a 
pseudogapped superfiuid, and of (b) T c as a function of 1 /kfU in a trap, where 
kf is the noninteracting Fermi momentum, and a is the s-wave scattering length. 



A detailed derivation of this result can be found in Ref. [34]. 
Note that the T-matrix is effectively a renormalized pairing in- 
teraction. It shares exactly the same pole structure as the two- 
particle Green's function, Gi. Through a Taylor expansion of 
its denominator, one can extract the pair dispersion: 



Q(Q)*z(ini-n q+t i pair ). 



(6) 



The superfiuid instability is given by 1 + Ux(0) = oc p pair , 
which is the BEC condition for pairs. Note that^(0 involves a 
mix of bare and full Green's functions. We emphasize that this 
is a natural consequence of the equation of motion technique 
since it involves the operator Gr . It is this GoG form ofx that 
leads back to the BCS-form of gap equation in the absence of 
finite momentum pairs. 

We focus on the superfiuid phase where t pg (Q) diverges at 
g = 0. Defining 



Q*0 



(7) 



we have 



G (-K) + cffi 



= - A z pg G (-K) + 5Z. 



(8) 



Y, pg takes the same form as T, s , 



Neglecting the residue term 8L, 
Thus we have immediately the BCS form of total self energy, 
Z(K) = -A 2 G (-K), with A 2 = A 2 C + A 2 pg . This then leads to 
the BCS form of gap equation, 



l-2/(£ k ) 
2E k 



= 0, 



(9) 



where = ~<J(ek - p) 2 + A 2 is the quasiparticle dispersion. 
Different from the BCS mean-field theory, we emphasize that 
here A 2 contains contributions from both condensed and non- 
condensed pairs so that it in general does not vanish at T c . 
Note that the finite q pairs are different from the order param- 
eter collective modes; the latter represent a coherent motion of 
the condensate. Here A 2 C and A 2 g are loosely proportional to 
the density of condensed and noncondensed pairs, respectively. 



Equations ^ and (0, along with the number equation 
» = 2^00*0, (10) 

K 

form a closed set of equations for the homogeneous case, which 
can be used to solve for p, T c , and the gaps at T < T c . T c is 
determined by setting A 5C = 0. Typical behaviors of the gaps 
are shown in Fig. 2(a). 

To address Fermi gases in a trap, we use the local density 
approximation, by replacing p. — » p. - V, rap (r). Then the num- 
ber equation becomes N - J d 3 rn(r). In Fig. 2(b) we show 
the BCS-BEC crossover behavior of T c in a trap. Here 1/kfa 
parametrizes the interaction strength. 

The RF response can be derived using the linear response 
theory. The RF interaction is described by 



H r f = e' n ' ^ d i x(j/^(j/ 1 + h.c, 
and the response Kernel by 

D(iCl,) = Yj G{2)(K)G(3)(K+Q) - 



(11) 



(12) 



We assume hyperfine level 3 was initially empty. In the ab- 
sence of final state interactions, as in 40 K, we obtain l35ll the 
RF current 

1 



I(k, v) = — Im D R {v + p-p 3 ) 



i-J]A(k, W )/M 



(13) 



In order to address A(k, of) = -2 Im G(k, to + /0 + ) properly, we 
need to include the lifetime effects of finite momentum pairs 
and add an incoherent term iLq in (and only in) J, pg , reflecting 
the residue term 51. which we drop in solving the set of equa- 
tions, i.e., 



E„„(k, a>) = 



pg 



a) + €k - p + iy 



- /En 



(14) 



While above T c the spectral function with a pseudogap consti- 
tute a double peak structure with suppressed spectral weight at 
the Fermi level, below T c , there is a zero at co — -(et - p). As 
A sc increases with decreasing T below T c , the spectral peaks 
sharpen rapidly. This is a phase coherence effect. The parame- 
ters y and Zo can be estimated from experimental RF spectra. 

4. Comparison between theory and experiment 

In Fig. 3, we compare between theory (curves) and experi- 
ment (symbols) (a) the density profile [36] and (b) system en- 
ergy IHD for 6 Li in the unitary limit. Both experimental and 
theoretical density profiles are very smooth, in good agreement 
with each other. Alternative theories predicts either a kink at the 
edge of the superfiuid core or nonmonotonic radial and temper- 
ature dependence. The energy comparison also reveals a quan- 
titative agreement. The fact that the unitary and noninteracting 
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Figure 3: Comparison of (a) density profile and (b) energy E/Ef for a unitary 
6 Li gas between theory (curves) and experiment (symbols). Also shown in (b) 



is comparison for the noninteracting energy. Here Ef 
acting Fermi energy. 



■ &h7Y is the noninter- 



curves merge at T « 0.67V » T c manifests the presence of a 
pseudogap. It should be noted that there is no fitting parameter 
in our theoretical calculations. 

Shown in Fig. 4 is a comparison of the spectral intensity map 
as a function of k and single-particle energy to + /u between 
experiment ll3Cill and theory j35ll for a unitary 4l) K gas at a tem- 
perature slightly above T c . The white dashed curve is the exper- 
imentally extracted quasiparticle dispersion, whereas the solid 
curve is obtained theoretically following the experimental pro- 
cedure. It is evident that theoretical and experimental results are 
rather close to each other. Indeed, as T decreases from above 
to below T c , the spectral intensity map evolves lL3~5ll from an 
upward dispersing branch at high T to a bifurcation around T c , 
and finally to a downward dispersing branch at T <k T c . This 
result establishes the actual single particle dispersions which 
contribute to the RF current, revealing that the broad peak in 
previous momentum-integrated RF spectra II 1 3fl indeed has a 
pairing origin. Furthermore, it also shows that, despite the trap 
inhomogeneity, momentum resolved RF spectroscopy can still 
provide a quantitative measure of the spectral function and sin- 
gle particle dispersion. It also lends support for the present GqG 
scheme since alternative NSR-based theories do not [37] seem 
to generate the two-branch-like feature observed in Ref. l30ll . 
The downward dispersion around (and above) T c provides di- 
rect evidence for the existence of a pseudogap above T c at uni- 
tarity. Our theory serves as a basis for momentum-resolved RF 
spectroscopy analysis. 

In summary, we have presented a pairing fluctuation theory 
where finite momentum pairing plays a progressively more im- 
portant role as the pairing strength increases, leading to a pseu- 
dogap in the single particle excitation spectrum. This theory 
has been successfully applied to multiple experiments in atomic 
Fermi gases. 

This work was supported by Zhejiang University and NSF of 
China Grant No. 10974173. 



References 

[1] Q. J. Chen, J. Stajic, S. N. Tan, and K. Levin, Phys. Rep. 412, 1 (2005). 
[2] J. N. Milstein, S. J. J. M. F. Kokkelmans, and M. J. Holland, Phys. Rev. 

A 66, 043604 (2002). 
[3] D. M. Eagles, Phys. Rev. 186, 456 (1969). 

[4] A. J. Leggett, in Modern Trends in the Theory of Condensed Matter 
(Springer- Verlag, Berlin, 1980), pp. 13-27. 



5 10 
k(um- 1 ) 



Figure 4: Comparison of spectral intensity map l(k, v)k 2 /(2n 2 ) between exper- 
iment (left) from Ref. 1 30] and theory (right). The white dashed curve is an 
experimental extracted quasiparticle dispersion, and the white solid line is ob- 
tained theoretically following the same experimental data analysis procedure. 



[5] P. Nozieres and S. Schmitt-Rink, J. Low Temp. Phys. 59, 195 (1985). 

[6] Y. J. Uemura, Physica C 282-287, 194 (1997). 

[7] R. Friedberg and T. D. Lee, Phys. Lett. A 138, 423 (1989). 

[8] Randeria, Physica B 198, 373 (1994). 

[9] B. Janko, J. Maly, and K. Levin, Phys. Rev. B 56, Rl 1407 (1997). 
[10] Q. J. Chen, I. Kosztin, B. Janko, and K. Levin, Phys. Rev. Lett. 81, 4708 
(1998). 

[11] C. A. Regal, M. Greiner, and D. S. Jin, Phys. Rev. Lett. 92, 040403 (2004). 

[12] Q. J. Chen, C. A. Regal, M. Greiner, D. S. Jin, and K. Levin, Phys. Rev. 
A73, 041601 (2006). 

[13] C. Chin, M. Bartenstein, A. Altmeyer, S. Riedl. S. Jochim, J. Hecker- 
Denschlag, and R. Grimm, Science 305, 1 128 (2004). 

[14] M. W. Zwierlein, C. A. Stan, C. H. Schunck, S. M. F. Raupach, A. J. 
Kerman, andW. Ketterle, Phys. Rev. Lett. 92, 120403 (2004). 

[15] J. Kinast, A. Turlapov, and J. E. Thomas, arXiv:cond-mat/0409283 

[16] Q. J. Chen, J. Stajic, and K. Levin, Phys. Rev. Lett.95, 260405 (2005). 

[17] J. Kinast, A. Turlapov, J. E. Thomas, Q. J. Chen, J. Stajic, and K. 
Levin, Science 307, 1296 (2005), published online 27 January 2005; 
doi:10.1126/science.l 109220. 

[18] M. W. Zwierlein, J. R. Abo-Shaeer, A. Schirotzek, and W. Ketterle, Na- 
ture 435, 170404 (2005). 

[19] Q. J. Chen, Y. He, C.-C. Chien, and K. Levin, Phys. Rev. A 74, 063603 

(2006) . 

[20] Y. He, C.-C. Chien, Q. J. Chen, and K. Levin, Phys. Rev. A75, 021602(R) 

(2007) . 

[21] P. Fulde and R. A. Ferrell, Phys. Rev. 135, A550 (1964); A. I. Larkin and 
Y. N. Ovchinnikov, Zh. Eksp. Teor. Fiz. 47, 1 136 (1964) [Sov. Phys. JETP 
20, 762 (1965)]. 

[22] G. B. Partridge, W. Li, R. I. Kamar, Y. A. Liao, and R. G. Hulet, Science 
311, 503 (2006). 

[23] M. W. Zwierlein, C. H. Schunck, A. Schirotzek, and W. Ketterle, Nature 

(London) 442, 54 (2006). 
[24] C. H. Schunck, Y. Shin, A. Schirotzek, M. W. Zwierlein, and W. Ketterle, 

Science 316, 867 (2007). 
[25] J. Kinnunen, M. Rodriguez, and P. Torma, Science 305, 1 131 (2004). 
[26] Y. He, Q. J. C hen, and K. Levin , Phys. Rev. A 72, 011602(R) (2005). 
[27] E. J. Mueller. larX"iv:07 1 1 .01 82 K unpublished). 
[28] S. Basu and E. Mueller, Phys. Rev. Lett. 101, 060405 (2008). 
[29] Y. He, C. C. Chien, Q. J. Chen, and K. Levin, Phys. Rev. Lett. 102, 020402 

(2009). 

[30] J. T. Stewart, J. P. Gaebler, and D. S. Jin, Nature (London) 454, 744 

(2008) . 

[31] G. B. Partridge, K. E. Strecker, R. I. Kamar, M.W. Jack, and R. G. Hulet, 

Phys. Rev. Lett. 95, 020404 (2005). 
[32] Q. J. Chen and K. Levin, Phys. Rev. Lett. 95, 260406 (2005). 
[33] C.-C. Chien, Q. J. Chen, Y. He, and K. Levin, Phys. Rev. Lett. 97, 090402 

(2006). 

[34] Q. J. Chen, Ph.D. thesis, University of Chicago, 2000, (available in the 

ProQuest Dissertations & Theses Database online). 
[35] Q. J. Chen and K. Levin, Phys. Rev. Lett. 102, 190402 (2009). 
[36] J. Stajic, Q. J. Chen, and K. Levin, Phys. Rev. Lett. 94, 060401 (2005). 
[37] E. Mueller, private communication. 



4 



